Quantum molecular dynamics simulations of thermophysical properties of fluid ethane 
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We have performed first-principles molecular-dynamics simulations based on density-functional 
theory to study the thermophysical properties of ethane under extreme conditions. We present new 
results for the equation of state of fluid ethane in the warm dense region. The optical conductivity is 
calculated via the Kubo-Greenwood formula from which the dc conductivity and optical reflectivity 
are derived. The close correlation between the nonmetal-metal transition of ethane and its decom- 
position, that ethane dissociates significantly into molecular and/or atomic hydrogen and some long 
alkane chains, has been systematically studied by analyzing the optical conductivity spectra, pair 
correlation functions, electronic density of states and charge density distribution of fluid ethane. 
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I. INTRODUCTION 



Pressure-induced transformations of hydrocarbons un- 
der extreme conditions have recently drawn exten- 
sive attention and have made giant achievement in 
astrophysics 1 . Specially, the planetary model with non- 
axisymmctric magnetic fields generated within a thin 
fluid shell has successfully revealed the non-axisymmetric 
nature of the Uranian and Neptunian magnetic fields^. 
As the interiors of these planets are composed of the 
fluids containing species, such as, C, H, N, O, etc, the 
knowledge of physical properties of fluid hydrocarbons 
(e.g., methane, ethane and benzene) under extreme con- 
ditions, such as equation of states (EOS), electrical con- 
ductivity and structure, is important in this context. 

Saturated hydrocarbon ethane (C2H6) has been re- 
ported as one of the most important decomposition prod- 
ucts of the giant planets^—. For instance, Richters et alA 
have reported that ethane is formed in the decomposition 
process of methane at high temperature (T=4000 K) and 
high pressure (P«100 GPa). The observation of atmo- 
spheric ethane in Neptune, where it might be brought up 
from the deep interior by convection process, could be ex- 
plained by the formation of ethane. Quantum molecular 
dynamics (QMD) simulations of Li et al^ also indicated 
that methane, which is a major constitute of the "ice" 
layer in Uranus and Neptune, mainly converts into ethane 
at density of 1.5 g/cm 3 along the principal Hugoniot. 
Meanwhile, ethane is an important fluid industrially, and 
the second member of the vitally interesting alkane series. 
Its chemical kinetics and thermophysical properties are of 
great interest for gas-turbine engines, high-speed propul- 
sion, and materials synthesis^. Pyrolysis and oxidation 
of ethane have been studied behind reflected shocked 
waves^. The EOS of ethane with temperatures ranging 
from 90 K to 625 K and pressures less than 70 MPa has 
been obtained by using Benedict- Webb- Rubin equation 9 . 
However, such studies on ethane are limited to narrow 
temperature and pressure region. More importantly, the 
nonmetal-metal transition (NMT) of fluid hydrocarbons 



is a crucial and hot topic recently. For instance, methane 
and benzene have been reported to transform to metal- 
lic fluids under extreme conditions, where the NMT are 
associated with the rapid C-H bond breaking and the for- 
mation of atomic or molecular hydrogen^^. The NMT 
of the system has also been found to significantly affect 
the optical properties, such as the reflectance and absorp- 
tion coefficient. The optical reflectivity of fluid deuterium 
has been successfully detected in dynamic compression 
experiments^. To date, the thermophysical properties of 
ethane under extreme conditions have not been reported 
yet. Considering the present facts that (i) ethane also de- 
notes a typical prototype and plays an important role in 
hydrocarbons, (ii) ethane is one of the most important 
decomposition products in the extreme-condition mat- 
ter evolution process of the giant planets, and (iii) as a 
same kind of alkane species, the differences and relations 
of the chemical dissociation properties between methane 
and ethane are interesting, it is thus highly necessary to 
study the equation of state, chemical decomposition and 
recomposition picture, electronic and optical properties 
of fluid ethane under extreme conditions by using an ef- 
ficient method. 

Materials under extreme conditions are complex, 
where partial dissociation, ionization and electronic de- 
generation coexist. QMD simulations, based on finite 
temperature density functional theory (FT-DFT)^ 2 -, of- 
fer a powerful tool to study thermophysical properties of 
complex materials under such extreme conditions^—. 
Combined with the Kubo-Greenwood formulatio n 17 ' 18 , 
the electrical and optical properties of various systems 
can be accurately calculate d 19 ' 20 . 

In the present work, we have performed QMD simu- 
lations to study the thermophysical properties of fluid 
ethane. By analyzing pair correlation functions, atomic 
structure, and charge density distribution, we have pre- 
dicted the chemical picture of the shocked ethane. In 
particular, the NMT for fluid ethane can be directly as- 
sociated to its dissociation, which is analogous to that of 
molecular fluid methane. Additionally, the optical prop- 
erties of fluid ethane are also discussed in the chemical 
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picture. 



II. COMPUTATIONAL METHODS 
A. First- principles molecular dynamics 

In this study, we have performed simulations for ethane 
by employing the Vienna ah initio simulation package 
(VASP) plane-wave pseudopotential code, which was de- 
veloped at the Technical University of Vienn a 2 1 i 22 . Elec- 
trons are fully quantum mechanically treated through 
solving the Kohn-Sham equations for a set of orbitals 
and energies within a plane- wave FT-DFT formulation^, 
where the electronic states are populated according to 
the Fermi-Dirac distribution at temperature T e . The 
ion-electron interactions are represented by a projec- 
tor augmented wave pseudopotential. The exchange- 
correlation functional is determined by generalized gradi- 
ent approximation with the parametrization of Perdew- 
Burke-Ernzerho^ 3 -. Atoms move classically according to 
the forces originated from the interactions of ions and 
electrons. The isokinetic ensemble (NVT) is employed by 
Nose-Hoover thermostat and the system is controlled 
in local thermodynamical equilibrium by setting the elec- 
tron temperature T e and the ion temperature Tj to be 
equal. Pseudopotentials with cutoff radius of 1.10 Bohr 
for carbon and hydrogen are adopted, and 650 eV plane- 
wave cutoff energy is employed in the simulations. 

In all the simulations, a total number of 128 atoms (16 
C2H6 molecules) are included in a supercell with periodic 
boundary condition. We perform finite-temperature, 
fixed- volume molecular dynamics simulations for selected 
densities ranging from 0.7 to 2.8 g/cm 3 along 3000, 4000 
and 8000 K isotherms. For molecular dynamic simula- 
tions, only r point of the Brillouin zone is included. All 
the dynamic simulations last 10 ps with time steps of 0.1- 
1.0 fs according to different conditions. For each pressure 
and temperature, the system is equilibrated within 1-2 
ps. The isotherm equation of states data and pair corre- 
lation functions are obtained by averaging over the final 
5 ps simulations for all the particles. 
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where ^k is Kohn-Sham eigenstate, with corresponding 
to eigenvalue e^.k, and occupation number /(e,-,k). tt is 
the volume of the supercell. W(k) is the K-point weight- 
ing factor. The i and j summations range over N discrete 
bands included in the calculation, and a over the three 
spatial directions. 

The imaginary part of conductivity 02 is obtained via 
the Kramers-Kronig relation as 
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where P is the principal value of the integral and v is 
frequency. The complex dielectric function e follows im- 
mediately from the two parts of the conductivity by the 
following equations 

ei(cj) = 1 a 2 (co); e 2 (co) = — cti(w). (3) 

LO LO 

The real part n(w) and the imaginary part k(co) of the 
refraction index have the following relation with the di- 
electric function 



e(co) — ei(io) + ie2(uj) — [n(u)) + ik(u)] 



(4) 



From these quantities, the absorption coefficient a(to) 
and reflectivity r(uS) are derived 



a(ui) = 
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III. RESULTS AND DISCUSSIONS 
A. Equation of state and structure 



B. Optical and electronic properties 

The electronic properties of dense ethane along 
isotherms are calculated in a spatial configuration of all 
atoms at a single time step within a MD trajectory. For 
dispelling the correlation time, the electronic conductiv- 
ity are evaluated for 15-20 selected snapshots taken from 
the QMD trajectory. The selected snapshots or config- 
urations are spaced by at least the correlation time. A 
4x4x4 Monkhorst-Pack k point mesh 2 ^ has been used 
to calculate the electronic structure. The dynamic con- 
ductivity is derived from the Kubo-Greenwood formula 



TABLE I: The simulated data of the density, pressure, and 
the standard deviation of pressure. 



p (g/cm 3 ) T (K) P (GPa) AP (GPa) 

0.80 3000 7.94 1.04 

0.90 3000 11.24 2.05 

1.00 3000 14.57 2.21 

1.20 3000 24.84 2.54 

1.40 3000 38.72 2.99 

1.60 3000 57.66 3.49 
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FIG. 1: (Color online) Pressure versus 
isotherms (T=3000, 4000 and 8000 K). 
of the simulation results along isotherms 
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The EOS has shown a systematic behavior in terms of 
density and temperature. As has been reported previ- 
ously, at temperatures between 3000 and 4000 K along 
the principal Hugoniot, the dissociation of methane oc- 
curs and nonmetal-metal phase transition takes place 
simultaneously 5 . As a consequence, we have studied the 
thermal EOS of fluid ethane along 3000, 4000 and 8000 K 



isotherms, as shown in Table I. The pressure P consists 
of two components P e and Pi. The gradient of total en- 
ergy calculated by VASP determines P e , which represents 
the contributions from ion-ion, ion-electron and electron- 
electron interactions. Due to the classical movement of 
ions, the ionic part Pi is derived by the ideal gas expres- 
sion. We thus have P = P e + p 71 KbT, where p n is the 
number density. The simulated isotherm EOS for dense 
ethane is shown in Fig. 1. We have fitted the simulated 
results for the pressure P by expansions in terms of den- 
sity p and temperature T as following P — ^A i jp l T^ . 
The coefficients Aij are listed in Table II. This expan- 
sion can easily be applied in hydrodynamic simulations 
for warm dense ethane in astrophysical applications. 

To clarify the structural transitions of ethane under ex- 
treme conditions, we have calculated the pair-correlation 
functions (PCFs), which represent the possibility of find- 
ing a particle at a distance r from a reference atom. The 
PCFs and the atomic structure together with charge den- 
sity distributions along 3000 K isotherm are presented in 
Fig. 2 and 3, respectively. At the low density p=1.60 
g/cm 3 , the peak in the C-H PCF gc-n(r) exists at 
around 1.09 A (Fig. 2a), which is the equilibrium in- 
ternuclear distance of the C-H bond in ethane molecule. 
Meanwhile, the peak in the C-C PCF gc-c( r ) occurs at 
the equilibrium internuclear distance of C-C bond of sat- 
urated hydrocarbon (1.54 A). Therefore, below the den- 
sity of 1.60 g/cm 3 , ethane remains in its ideal molecular 
configuration without dissociation. From the structure 
of ethane at 1.60 g/cm 3 (Fig. 3a), we can straightly find 
that all of the ethane molecules are not dissociated. At 
the density of 2.00 g/cm 3 , we find a significant reduc- 
tion of the maximum of gc-nif) around 1.09 A, which 
indicates that the C-H bond breaks rapidly and ethane 
molecules dissociate with the increase of the density. On 
the other hand, around 0.75 A a shoulder appears in the 
H-H PCF at this density, which implies the formation of 
molecular hydrogen. As the density further increases to 
2.3 g/cm 3 , the maximum of the PCF <7c-h(?") continues 
to reduce while the peak amplitudes in <?c-c( r ) increases, 
which indicates that ethane further dissociates into some 
hydrocarbons. From Fig. 3c, we can see that a new 
species of long alkane chains show up in the atomic con- 
figuration, corresponding to the increase of the peak of 
.9c-c( r )- More hydrogen molecules forms at this density, 
which corresponds to the larger value of <7h-h(?") around 
0.75 A. Note that hydrogen molecules or atoms appear 
much more in Fig. 3c than in Fig. 3b. However, the 
peak of gH_-n(r) around 0.75 A gets less apparent, and 
even vanishes at 2.8 g/cm 3 . This may be because that 
(i) the peak of <7h-h(7) around 0.75 A are covered by the 
widening of the peak around the distance between the 
two hydrogen atoms bonding with a carbon atom (i.e., 
more complicated alkane structures and more hydrogen 
molecules and atoms form); and/or (ii) the formed hydro- 
gen molecules partially dissociated into hydrogen atoms. 

From the reported of Li et al£, it is known that above 
the density 1.05 g/cm 3 along the Hugoniot curve (above 
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TABLE II: Pressure (GPa) expansion coefficients Aij in terms 
of density (g/cm 3 ) and temperature (K). 



27.3780 -73.3830 53.4723 

1 0.0017 -0.0018 0.0020 
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FIG. 3: (Color online) The atomic structure and charge den- 
sity distribution of dense ethane at four density along 3000 
K isotherm, carbon and hydrogen atoms are denoted by red 
and blue balls, respectively, (a) p=1.60 g/cm 3 ; (b) p=2.00 
g/cm 3 ; (c) p=2.30 g/cm 3 ; (d) p=2.80 g/cm 3 . 



B. Optical and electronic properties 



FIG. 2: (Color online) Pair-correlation functions for H-H 
(black solid line), C-C (red dashed line), and C-H (blue dashed 
dot line) for four densities of dense ethane along 3000 K 
isotherm, (a) p=l.W g/cm 3 ; (b) p=2.00 g/cm 3 ; (c) p=2.30 
g/cm 3 ; (d) p=2.80 g/cm 3 . 



the pressure 25 GPa), methane begins dissociate and 
leads to the appearance of ethane and some linear carbon 
chains. In our study on the PCF of ethane, we find that 
ethane begins to dissociate at the pressure 81 GPa, and 
hydrogen and alkane chains form. Compared with the 
dissociation of methane at warm dense region^, both of 
the hydrocarbons have the similar dissociation behavior 
with increase of the density. The dissociated products 
of both alkane species are hydrogen and alkane chains. 
Therefore, starting at methane or ethane should result 
in basically the same final state comprising alkane chains 
and hydrogen molecules and atoms at high enough pres- 
sures and temperatures, except for the quantity of alkane 
chains and molecular or atomic hydrogen. In addition, 
in the whole dissociation process of ethane and methane 
simulated theoretically, no diamond is formed, whereas it 
was experimentally reported that the diamond is formed 
at pressure P=2Q GPa and temperature T=2000 K for 
methane 1 . Richters et alA have employed a homogenous 
nucleation mechanism, which is similar to the proposal 
of the direct graphite-to-diamond transitio n 26 ' 27 , to ex- 
plain the discrepancy between theory and experiment for 
methane. 



We have calculated the electronic conductivity and cor- 
responding reflectivity of fluid ethane along isotherms 
by Kubo-Greenwood formula. The behavior of the 
frequency-dependent conductivity u\{ui) at different den- 
sities along 3000 K isotherms is shown in Fig. 4. It 
has been found that o\ (w) at these four different densi- 
ties exhibit a uniform feature that their peaks all locate 
around 10.0 eV, which is correlated with the transitions 
to the lowest excited states. With the increase of den- 
sity along isotherm, the main peak of the electrical con- 
ductivity o\{lo) increases in amplitude and moves to low 
frequency. Such change of the peak indicates another 
quantity, dc conductivity, given as <7d c = limci (cj), in- 

creases with density. The dc conductivity is extracted 
from dynamic conductivity and shown in the inset of Fig. 
4. The dc conductivity a^c at low density along 3000 K 
isotherm is small, which corresponds to the character of 
a semiconducto r 28 ' 29 . With increasing density, the dc 
conductivity er^c increases rapidly. Above 2.30 g/cm 3 
(P=165 GPa), a metal-like conductivity is found. In or- 
der to analyze the effect of temperature on the nonmetal- 
metal transition, we have also shown the dc conductivity 
along higher temperature 4000 and 8000 K isotherms in 
Fig. 4. At low density and high temperature the dc con- 
ductivity can also reach a large value, which is a typical 
character of metallic behavior. 

Mott et al£2- have reported a decisive physical asser- 
tion, which is used to define the metallicity of a disor- 
dered system. When the characteristic mean free path 
of the valence electrons exceeds the mean distance be- 
tween the constituent atoms or molecules, the states in 
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FIG. 4: (Color online) Optical conductivity spectra along the 
3000 K isotherm. Data have been averaged over 15-20 uncor- 
rected MD configuration. The dc conductivities as a function 
of density along 3000, 4000 and 8000 K isotherms have been 
shown in inset. 



FIG. 6: (Color online) Variation of the frequency-dependent 
reflectivity r(oj) along 3000 K isotherm, variation of the re- 
flectivity as a function of density at fixed wavelengths of 404 
nm along 3000, 4000 and 8000 K isotherms. 
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FIG. 5: (Color online) The density of state of fluid ethane 
at four different density points along 3000 K isotherm. The 
Fermi energy is set to be zero. 



disordered system with any high-temperature will remain 
metallic or indeed attain metallic status. The constituent 
atoms or molecules provide the carriers of the electrical 
current. Applying Mott's criterion for the minimum elec- 
trical conductivity for fluid metals, the minimum elec- 
trical conductivity is about 1000 fi _1 cm _1 for fluid hy- 
drogen, rubidium, caesium and mercury 3 ^. The NMT 
has recently been observed in the nominally nonmetallic 
chemical elements hydrogen, nitrogen, and oxygen under 
extreme condition o 31 ' 32 . 

At high pressures, materials undergo dissociation 
bringing with some new interesting properties. For in- 
stance, the band gap between the valence and conduction 
bands decreases with pressure and then disappears. The 
NMT has close relation with the molecular dissociation 
in the shocked fluid. In the case of fluid hydrogen and 
methane, QMD simulations have suggested that high- 
pressure NMT and the dissociation of molecules with in- 
creasing pressure are closely connected^ 3 ^. The dissoci- 



ated atoms act as dopants and progressively occupy the 
dense-fluid band gap. In order to validate this physical 
picture, we have calculated the various electron density 
of states at different densities along 3000 and 4000 K 
isotherms as shown in Fig. 5. For the low density 1.60 
g/cm 3 , a gap exits, locating near the Fermi energy where 
only thermally activated electron transport occurs, which 
is exactly the situation in semiconductors or nonmetals. 
With increasing the density, the gap is gradually reduced, 
and finally disappears. At density 2.30 g/cm 3 , the elec- 
tronic states fill up the gap so that metal-like conductiv- 
ity follows. In addition, Compared with the same density 
along different isotherm, we can see the temperature ef- 
fect on the electronic structure. Such effect can also be 
visualized in the dc conductivity ddc-pressure plotted in 
the inset of Fig. 4. 

The optical reflectivity versus photon energy at dif- 
ferent densities along the 3000 K isotherm is shown in 
Fig. 6 (the main panel). With the increase of density 
along the isotherm, the reflectivity asymptotically ap- 
proaches unity for zero photon energy. The reflectivity 
versus density at fixed wavelength 404 nm (correspond- 
ing to 3.075 eV) of fluid ethane along the 3000, 4000, and 
8000 K isotherms have been shown in the inset of Fig. 6. 
As is seen, a measurable reflectivity arises from 0.10 to 
0.49-0.59 at typical wavelength of 404 nm spanning the 
visible spectrum, which is related with the high-pressure 
NMT. Similar behavior in the optical reflectance of fluid 
methane has been obtained by using QMD simulations 5 . 
The results of the optical reflectivity can be inspected 
through future experiments, and the reflectivity of liquid 
deuterium has been executed in experiment^. 



IV. CONCLUSIONS 

In summary, first-principles molecular-dynamics simu- 
lations based on density-functional theory have been used 
to study the thermal EOS and nonmetal-metal transition 
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of fluid ethane under extreme conditions. Systematic 
descriptions of pair-correlation function, atomic struc- 
ture, and the charge density distribution with various 
densities along the isotherms are obtained to show the 
dissociation process of molecular fluid ethane, which is 
associated with the nonmctal-metal transition. Ethane 
dissociate with increasing the density and temperature 
into molecular and/or atomic hydrogen and some long 
alkane chains, which is similar to the case of methane. 
However, no diamond-like carbon microstructure is ob- 
served for ethane under extreme condition, the present 
result shows a different high pressure behavior compared 



with previous experiments. The nonmetal-metal tran- 
sition has significant influence on the optical properties 
of fluid ethane, which needs to be verified in dynamic 
compression experiments. Our results are expected to be 
revealing for the planetary models of Neptune. 
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